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Abstract. We apply the method of nonUnear steepest descent to compute 
the long-time asymptotics of the Korteweg— de Vries equation with steplike 
initial data. 



1. Introduction 

We study the long-time asymptotic behavior of solutions of the Korteweg-de 
Vries (KdV) equation 

(1.1) qt{x,t) = 6q{x,t)qx{x,t) - qx.j;x{x,t), {x,t) eRxR, 

with steplike initial data q{x,0) = q{x) E C^^{R) such that 

q{x) — 0, as X — ?► +oo, 

q{x) — >■ — c^, as a; — >■ — oo, 



(1.2) 
moreover, 

(1.3) / e'^"''{\q{x)\ + \q{-x) + c'^\dx<oo, Co > c> 0, 



(1.4) {x^ + l)\q^'^x)\dx <oo, i = l,...,n. 
Jr 

It is known (cf. |15| . |16|). that this Cauchy problem has a unique solution satisfying 
g(-,t) e C^iR) and 

p+oo 

(1.5) / \x\{\qix,t)\ + \q{~x,t) + c^\)dx <oo, teR. 
Ja 

From several results ([3]-[8], [19], [20], [17], [29]), obtained on a physical level of 
rigor, it is known that the asymptotic behavior of q{x, t) as i — >■ oo can be split into 
three main regions: 

• In the region x < —6c^t the solution is asymptotically close to the back- 
ground — up to an decaying dispersive tail. 

• In the region — 6c^i < x < Ac^t the solution can asymptotically be de- 
scribed by an elliptic wave. 

• In the region < 4c^t < x the solution is asymptotically given by a sum of 
solitons. 



2000 Mathematics Subject Classification. Primary 37K40, 35Q53; Secondary 37K45, 35Q15. 

Key words and phrases. Riemann— Hilbert problem, KdV equation, steplike. 

Research conducted in the framework of the project "Ukrainian branch of the French- Rus- 
sian Poncelet laboratory" — " Probability problems on groups and spectral theory" . Research 
supported by the Austrian Science Fund (FWF) under Grant No. Y330. 

1 



2 



I. EGOROVA, Z. GLADKA, V. KOTLYAROV, AND G. TESCHL 




Figure 1. Numerically computed solution q{x,t) of the KdV 
equation at time t = 10, with initial condition q{x, 0) = ^{eri{x) — 
1) — 5sech(a; — 1). 



This is illustrated in Figure [TJ In fact, the long-time asymptotics for this problem 
were first studied by Gurevich and Pitaevskii [19], [20]. These authors have used 
the Whitham multi-phase averaging method and obtained the main term of the 
asymptotics of the solution in terms the Jacobi elliptic function. Moreover, they 
gave a qualitative picture of the splitting of an initial step into solitons. Since the 
Schrodinger operator with the Heaviside step function as potential has no discrete 
spectrum, this picture refuted the general idea that solitons arise only from the 
discrete spectrum. This phenomenon was explained by Khruslov j21| . |22| with 
the help of the inverse scattering transform (1ST) in the form of the Marchenko 
equation. The 1ST not only made it possible to obtain an explicit form of these 
asymptotic solitons but also to give a rigorous proof that the solitons are generated 
by a small vicinity of the edge of the continuous spectrum. Further developments 
of this method can be found in [31] and [23] . The first finite-gap description of 
the asymptotics for the steplike initial problem of the KdV equation was given 
by Bikbaev and Novokshenov [6] only in 1987 (see also [3]-[5], [8], [7], and the 
review [29]). The results are based on an analysis of the Whitham equations and 
the theory of analytic functions on a hyperelliptic surface. Our aim here is to use 
the nonlinear steepest decent method for oscillatory Riemann-Hilbert problem (see 
[l8] for an introduction to this method in the case c = as well as for further 
references) and apply it to rigorously establish the above mentioned asymptotics. 
Related results for an expansive step (— > 0) can be found in [27] . 

The paper is organized as follows: Section [5] provides some necessary information 
about the inverse scattering transform on steplike backgrounds. Then we establish 
the asymptotics in the soliton region in Section |3j In Section |4] the initial RH 
problem is reduced to a "model" problem in the domain —6c^t < x < Ac^t, and in 
Section|5]we solve this model problem. Section |6| contains the solution of the model 
problem in the domain x < —Q<?t. 
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2. Statement of the RH problem and the first conjugation step 

Let q{x,t) be the solution of the Cauchy problem ()1.1|) - (|1.4|) . Associated with 
q{x,t) is a self-adjoint Schrodinger operator 

(2.1) H{t)^^^ + q{.,t), ^{H) ^ H\R) C L\R). 

Here L'^(R) denotes the Hilbert space of square integrable (complex- valued) func- 
tions over R and H''{M.) the corresponding Sobolev spaces. 

The spectrum of H consists of an absolutely continuous part [— c^, oo) plus a finite 
number of eigenvalues —Kj G (— oo, — c^), ^ < j < N , where c < ki < ... < kn. In 
turn, the absolutely continuous part of the spectrum consists of the part [0,cxd) of 
multiplicity two and the part [— c^,0] of multiplicity one. In addition, there exist 
two Jost solutions (l){k,x,t) and (j)i{k,x,t) which solve the differential equation 

(2.2) H{t)^{k, X, t) = k'^^ik, X, t), Im(fc) > 0, 

and asymptotically look like the free solutions of the background equations 

(2.3) lim e~"''^0(fc,a;,t) = 1, lim e''=i^(/)i(fc, x, t) = 1. 

a;— f-t-oo a;— ^ — oo 



Here ki — V^M-c?, and fci > for k e [0,ic)r. The last notation means the right 
side of the cut along the interval [0,ic]. Accordingly, ki < for k G [0,ic)i, i.e. 
from the left. As a function of k the function (j){k, x, t) (resp., 0i) is analytic in the 
domain ~ {k : Im(fc) > 0} (resp. := C*^ \ (0, ic]) and continuous up to the 
boundary of this domain. Here subscript U corresponds to the upper half plane. 

The Jost solutions admit the usual representation via the transformation oper- 
ators 

0(fc, X, t) = e"=^-H / K{x, y, ty'^vdy, 
(2.4) ■'^ 

</)i(fc,x,i)=e-'^-i^+ / Xi(x,y,i)e-"=^^dy, 



where K(x,y,t) and Ki{x,y,t) are real valued functions, and 
(2.5) K{x,x,t)^-J q{y,t)dy, K^{x,x,t) ^ - j {q{y,t)+^)dy. 
Furthermore, one has the scattering relations 



T{k,t)(j)i{k,x,t) ^(t){k,x,t) + R{k,t)(j){k,x,t), fc G M, 

(2.6) 

Ti{k,t)(j){k,x,t) ^(f)i{k,x,t) + Ri{k,t)(t)i{k,x,t), ki € M, 

where T(fc,t), R{k,t) (resp. Ti{k,t), Ri{k,t)) are the right (resp. left) transmission 
and reflection coefficients. They constitute the entries of the scattering matrix. 
Denote by 

(2.7) W{k,t) = (l)i{k,x,t)<p'{k,x,t) - (l)[{k,x,t)<p{k,x,t) 

the Wronskian of the Jost solutions, where /' = In what follows we assume 

that the initial data (|1.2p belong to the generic class of nonresonant potentials for 
which 



(2.8) 



W{±ic, 0) 7^ 0. 
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Lemma 2.1 ([ID], [E]). The entries of the scattering matrix have the following 
propertied: 

1. The transmission coefficients T{k,t), Ti{k,t) are meromorphic in the domain 

:— \ {0,ic], continuous up to the boundary and have simple poles at 
iKi, . . . ,iKN- The residues ofT{k,t) are given by 

(2.9) Resi,^. T{k,t) = {t)jj {t)\ where -fj{t)-^ = ||0(iK,, t)|i2, 

and (j){iHij,x,t) = fij(t)(j)i{iKj,x,t). 

2. Everywhere in the domain 

(2.10) T{k,t) =2ikW~^{k,t), Ti{k,t) ^2ikiW-^{k,t). 



3. The reflection coefficient has the symmetry property R{~k, t) — R{k, t) as k Cz 
and 



(2.11) Ti{k,t)T{k,t) = l-\R{k,t)\^, R{k,t)T{k,t) + Ri{k,t)T{k,t) = 0, fc e R; 

4. The functions i?i,Ti and T also possess the symmetry property with respect to 
ki £ R, in particular, Ti{k{—ki),t) = Ti{k{ki),t) . Moreover, 



(2.12) --Tik,t)T-^k,t) = Ti{k,t)T^\k,t) = Ri{k,t), kie[-c,c]. 

5. The time evolutions of the quantities Jj{t), R(k,t) and \T{k,t)\'^ are given by 
R{k,t) = i?(fc)e8''='* for k e R, \Ti{k,t)\'^ = |ri(A:)|2e8"='* for k e hic,ic], and 
jj{t) = 7je''''?*, where 7^ = 7^(0), R{k) = R{k,0) and T{k) = T{k,0). 

6. Under the assumption (|1.3p the function R{k) admits an analytic continuation to 
the domain D := {fc : | Inifc| < Co}\ [— ic, ic] preserving the symmetry properties 
R{~k) = R{k) for ken and i?(fc(-fci)) = i?(fc(fci)) for fci € [-c,c]. 

The properties, cited in this lemma, belong to the list of necessary and sufficient 
properties of the scattering data for the step-like potential with prescribed behavior 
of the perturbations. All of them, except of the last one, are valid for much wider 
class of perturbations, then the class (|1.3p . for example, for the class of potentials 
with the first moment of perturbations finite. 

Lemma 2.2. Abbreviate 



(2.13) xik)-=- hm^r(fc + £,0)Ti(fc + £,0) 
for k G ic, ic] . Then 

(2.14) x{-k) = -xik) and > for fce[-ic,ic]. 

Proof. Immediate from (|2.10l) . □ 

Consider a vector-function m(fc, x, t) as a function of spectral parameter k, k £ 
C\ [— ic, ic], where x,t are fixed parameters. We define this vector-function as 
follows 



(T(fc,i)</.i(fc,a;,t)e"=^, </>(fc, x, i)e-"=^) , fc G C^, 



(2.15) m{k,x,t)^^ (0(-fc,x,i)eK ^ T{-k,t)M-ix,t>~"'h , k € 
where C^:={fc: Im fc > 0} \ (0, ic], {fc : Im fc < 0} \ (0, -ic]. 



"'^We list here only those of the properties, that are relevant for the present paper 
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Lemma 2.3. The function m{k) = m{k,x,t), defined by formula ()2.15p has the 
following asymptotical behavior 



We are interested in the jump condition of m{k,x,t) on the contours S U Ec, 
where S = M, oriented LTR (left-to-right), and Sc = [ic, — ic], oriented top-down. 

In general, for an oriented contour S, the value m+(fc) (resp. m^{k)) will denote 
the nontangential limit of m{K) as k — fc from the positive (resp. negative) side of E. 
Here the positive (resp. negative) side is the one which lies to the left (resp. right) as 
one traverses the contour in the direction of its orientation. In order to not mix up 
limit values of functions from the different sides of contours with another meaning 
of signs -|- and — , in what follows we denote the upper (resp. lower) half plane as 
(resp. C^). Any notation, which is connected with upper or lower half plane, 
will be also marked by subscript U or L. For example, = [ic, 0]. Moreover, 
by subscripts I and r we will mark, when necessary, the values of functions from 
the left and right of the cut [— ic, ic]. For example, as the definition (see (j2.13p ) 
of the function x we could write x — or x = —[TTi]r. Note also, that 

the reflection coefficient i?(fc), k £ M., the function x{k), k G [— ic, ic] and the 
discrete spectrum together with right normalizing constants (Kj,7j), 1 < j < -/V, 
completely define the kernel of the right Marchenko equation, and, therefore, the 
potential q{x) (cf.ITU], [H])- That is why we refer to them as the minimal scattering 
data of operator H{0). 

Theorem 2.4. Let {R{k), k £ K;x(fc), k £ [-ic,ic]r; (Kj-,7j), I < j < N} be the 
minimal scattering data of the operator H{0). Then m{k) — m{k,x,t) defined in 
()2.15|) is a solution of the following vector Riemann-Hilbert problem. 

Find a vector-valued function m{k) which is meromorphic away from EUEc with 
simple poles at iiKj and satisfies: 

(i) The jump condition mjf.{k) — m^{k)v{k) 



(2.16) 




Proof. Will be given in Appendix \K\ 



□ 



(2.17) 




(ii) the pole conditions 



(2.18) 
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(Hi) the symmetry condition 



(2.19) 



m{—k) 



(iv) the normalization condition 



(2.20) lim miin) = (1 1). 



Here the phase $(A;) — ^{k,x,t) is given by 
(2.21) $(/s) = 8ifc^ + 2iA:|, 

Proof. Will be given in Appendix 1X1 



□ 



We note that m{z) defined in (|2.15l) is the only solution of the above Riemann- 
Hilbert problem. This can by seen after rewriting the pole conditions as jump 
conditions (see below) from 18, Thm. 3.2] (or alternatively from [28l Thm. 4.3]). 
Since all conjugation and deformation steps applied below are reversible, the solu- 
tions of all further Riemann-Hilbert problems will be unique as well. In this respect 
note that our jump matrix satisfies 



and det(w(fc)) = 1. 

For our further analysis we rewrite the pole condition as a jump condition 
and hence turn our meromorphic Riemann-Hilbert problem into a holomorphic 
Riemann-Hilbert problem following literally T^. Choose e > so small that the 
discs \k — \Kj\ < e lie inside the the domain and do not intersect any of the 
other contours. Denote the circle boundaries of these small discs as T^. Redefine 
m(fc) in a neighborhood of iKj respectively —iHj according to 



Note that in we redefined m{k) such that it respects our symmetry (|2.19l) . Then 
a straightforward calculation using ReSiKm(fc) = limfc^iK(fc — m)m{k) shows the 
following well-known result: 

Lemma 2.5 ( 18 ). Suppose m{k) is redefined as in ()2.23p . Then m{k) is holo- 



morphic m C \ (S U Sc U U^i(Tf U Tf)) . Furthermore it satisfies ([2T7l) . ([2l^ . 



(2.22) 




(2.23) 




(P:^ and 



(2.24) 
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where the small circle around iKj is oriented counterclockwise and the one around 
—iKj is oriented clockwise. 



3. ASYMPTOTICS IN THE DOMAIN 4c^i < X 

To reduce our RH problem to a model problem, that can be solved explicitly, we 
will use the well-known conjugation and deformation techniques. 

Lemma 3.1 (Conjugation). Let m be the solution of the RH problem m+(k) = 
m-{k)v[k), k £ H. Assume that S C E. Let D be a matrix of the form 



(3.1) Dik) 



(d{k)-^ 
d{k) 



where d : C\S — > C is a sectionally analytic function. Set 

(3.2) m{k) = m{k)D{k), 
then the jump matrix transforms according to 

(3.3) v{k) = D^{k)-^v{k)D+{k). 

If d satisfies d{k) ^ 0, d(~k) = d{k)^^ for € C \ E and limK-s-oo '^(i'^) = 1; then 
the transformation m{k) = m{k)D{k) respects the symmetry and normalization 
conditions (|2.19p and (|2.19p . respectively. 

In particular, we obtain 
respectively 

(^^^11 vi2d+d^\ ~ 

Now we make the first conjugation step, which allows us to take into account the 
influence of the discrete spectrum. To this end we will need the value kq defined 
via Re($(iKo)) = 0, that is, 

(3.6) ^„ = > 0. 

We will set kq = if f < for notational convenience. Then we have Re(<l'(iKj)) > 
for all nj > kq and Re($(iKj)) < for all kj > kq. Hence, in the first case the 
off-diagonal entries of our jump matrices are exponentially growing and we need to 
turn them into exponentially decaying ones. Therefore we set 



(3.7) A(fc):= n 



k 
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and introduce the matrix 
(3.8) 

1 



k—iK-j 



D{k) = < 



. 2 t<^>(ift^) 





k-\-iKj 



' I Do{k), \k-iK,\<e, j = l,...,N, 





'^l'"^ |i^o(A:), \k + iK,\<e, j = l,...,iV, 



else, 



where 



i?o(fc) = 

Observe that by A(— fc) = A^^(fc) wc have 



/A(fc)-i 
V Hk) J ' 



(3.9) 

Now we set 
(3.10) 



D{^k) 



1 

1 



D{k) 



1 

1 



m(fc) = m{k)D{k). 



Note that by p.9p this conjugation preserve properties (|2.19p and (|2.20p . 

Then (for details see Lemma 4.2 of [18]) the jump corresponding to kq < nj is 
given by 



(3.11) 



v{k) = 
v{k) = 



^0 



(fc-iKj)A^(/c) ' 



1 



fc e T' 



(7 



\ i7|c**<'"j'A2(fc) y 

and the jumps corresponding to kq > Kj (if any) by 



(3.12) 



v{k) 
w(/fc) 



1 



i7jc ^ 



„2„t*(i»,)A(fc)- 



i^2^'*(i'=j)A(fe) 

1 



.0 



fc G TV 



fc e Tf. 



In particular, all jumps corresponding to poles, except for possibly one if kj — kq, 
are exponentially close to the identity for t — > oo. In the latter case we will keep 
the pole condition for kj = kq which now reads 



(3.13) 



ReSiKj. TO(fc) — ^lim m(fc) | .^2 



Res_iKj TO(fc) = lim m(fc) 





-i72e**('''^)A(iKj)-2 
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Furthermore, the jump along E U Sc now reads 



(3.14) w(fc) = < 



l-|i?(fc)|2 -A2(/c)i?(fc)e-**('=) 
A-2(fc)i?(fc)e**W 1 



1 

A-2(fc)x(fc)e**« 1 

1 A2(A:)x(fc)c-**('^) 
1 

The new Riemann-Hilbert problem 

(3.15) m+(fc) = m_(fc)u(fc) 



u 



fc e s 



k e E^ 



for the vector m preserve its asymptotics (|2.20p as well as the symmetry condition 
(|2.19p . In particular, after conjugation all jumps corresponding to poles are now 
exponentially close to the identity as f — >■ oo. To turn the remaining jumps along 
S U Sc into this form as well we chose two contours S'^ U sufficiently close to 
the original contour such that R is analytic in the region ft^ U il^ between these 
contours (cf. Figure [J). 



IKl 




Figure 2. Contour deformation in the soliton region. 



Now we factorize the jump matrix along S according to 
(3.16) V = bl^bu 

and set 



1 -A\k)R{k)\ ( 1 
1 ) U-2(fc)i?(/c) 1 



(3.17) 



m{k)b^\k), keOF, 
m{k) = { m{k)bl\k), k G VL^, 
rh{k), else, 



such that the jump along S is moved to T,^ U and given by 



(3.18) 



v{k) = < 



1 







A-'^{k)R{k) 1 



keY.'- 



1 -A2(fc)i?(fc)\ ^^^^ 



10 I. EGOROVA, Z. GLADKA, V. KOTLYAROV, AND G. TESCHL 

The jumps along the circles U Tj are unchanged and the jump along Ec now 
reads 



(3.19) v{k) 



- i?+ + x)A-2 i)' "^^c, 



1 (i?_ -R+ + x)A^ 
1 



k e 



The following lemma shows that this jump in fact also disappears. 
Lemma 3.2. The following identities are valid: 

R^{k)- R+{k) + x{k) =0, fceE^, 

i?_(fc)-i?+(fc)+x(fc) =0, fceE^ 
Proof. With the help of the Pliicker identity (cf. [SO])- □ 

Hence, all jumps v are exponentially close to the identity as t — > oo and one can 
use Theorem A. 6 from 25 to obtain (follow literally the proof of Theorem 4.4 in 

Theorem 3.3. Assume p.3p - p.4p and abbreviate by Cj — 4k| the velocity of the 
j 'th soliton determined by Re($(iKj)) = 0. Then the asymptotics in the soliton 
region, x/t ~ Ac? > C for some C > Q, are as follows: 

Let e > be sufficiently small such that the intervals [cj — s, Cj + e], 1 < j < N , 
are disjoint and lie inside (4c^,oo). 

If \ j — Cj\ < e for some j, one has 



-AKj-f^{x,t) 

(l + (2^,)-i7^(x,t))2 



(3-20) qix, t) - , ,,,, + Oit-') 

for any I G N, where 



l^i "t" l^j 



N 

(3.21) ^l^^^^) ^^2g-2.,.-+8«ft -Q 

If \ j — Cj\ > e, for all j, one has 

(3.22) q{x,t) ^ 0{t-^) 
for any I G N. 

4. Reduction to the model problem in the domain -6c'^t < x < Ac^t 

Now we turn to the elliptic region —Qc^t < x < Ac^t we first proceed as in the 
previous section to obtain v where we now use 

(4.1) A{k) :^ n 

since clearly kj > c > kq for all j. For the conjugation step we will use a g-function 
as first outlined in [13] . Our approach here is similar to [5] and [21] ■ 
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Set C = ift, then $(fc) = «>(fc,^) = Sifc^ + 2Aik^. Following ^ in the domain 
! \ Sc introduce the function 



(4.2) 



g{k) ;= g{k, x, t) ^ 12 / (fc^ + n'^) 



u2 , ..2,J>^±^ 

k"^ + c? 



dk 



where the parameters a = a(i^), < a < c and /i = /i(^), < /i < a < c, 



(4.3) 

are chosen to satisfy conditions 



(4.4) 



dk^O 



and 

(4.5) 5(fc) - 4fc3 - ^ 0, k^oo. 

As is shown in these conditions can be satisfied for all values of parameter ^ 
in the domain — ^ < ^ < y . Set = [ia, — ia] with the orientation top-down. 

Lemma 4.1 ( [24j ) . The function g{k) possess the following properties 

(a) Function g is an odd function in the domain C\ Sc; g{k) — ^g{~k); 

(b) g_(fc) + g+{k) = as fc G \ E,,; 

(c) g_(fc) - g+(fc) = B as k eT,a, where B := B(C) = 2.g+(ia) > 0; 

(d) the asymptotical behavior holds as fc —> oo; 

(4.6) '-m - i,(fc, - ^^^^^^ - + ' ^"^^^^^^ 

2 2Ki 



Proof. The last property follows immediately from (I4.2p - (|4.5p and property (a). 

□ 

The signature table for the function Im g{k) is depicted in Figured 





— 


' IC 








k ia 




+ 






+ 














r^ia 






+ 


, — ic 





Figures. Signoflm(g) 



Introduce the function 
(4.7) d{k,t) =exp(t$(fc)/2-it5(fc)). 
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According to (|4.6|) we have 
(4.8) 



dik.t) = 1 + +0 
k\ 



fc3 



^(0 = 



12e(c2 - a(02) + 3c4 + 9a(0^ - 6a(0'c2 



Since the functions $(fc) and ^(/c) are both odd functions of fc, the function (i(-,t) 
is analytic in C \ Ec and satisfies d{—k,t) = d~^{k,t) plus d(k,t) — > 1 as fc — !■ oo. 
Let m(fc) be the solution of the problem p.lip - (l3.15p . Set rh{k) — fh{k)D{k,t), 
where the diagonal matrix D{k,t) is defined by (|3.ip with d{k,t), defined by (|4.7I) . 
Applying Lemma 13.11 we arrive at the following Riemann-Hilbert problem: 

rh^(k) — m^{k)v{k), m(fc) — )• (1, 1), fc — )• oo, 



(4.9) 
where 



(4.10) v{k) 



1 hnk,tt) 



hf{k,^,t) 1 



l-\R{k)\^ 
A-2(fc)i?(fc)e2«9W 

A-2(fc);^(fc)ei*(9++9- 



j = l,...,iV, 



,iV, 



-A2(fc)i?(fc)c-2«s(fc)\ 







ei*(s+-9-) A2(fc)x(fc)e-'*(9++f-' 
e~'*^^+~^~^ 



fees, 



where the entries 

(4.11) h^{k,tt) = i7r2(fc_iK^.)A2(fc)e*(*W-*('«^)-2i3(fc))^ 



(4.12) /ij^(fc,C,0 =i77'(fc + i'«j)A-2(fc)e-*(*W-*(-'«^)-2'9(fc))^ 

of the conjugation matrices on the circles decay exponentially with respect to t. 

Introduce two domains Vl^ and 17^, bounded by E and contours YP and E^ 
respectively, where the contours E*^ and E^ are symmetric with respect to map 
k —k and oriented LTR (cf. Figure |4]). Moreover, V.^ and fl^ must remain in 
the region where Im^g) > and Im(g) < 0, respectively. 

Following the standard procedure (see, for example, [H], [15], [53]) we factorize 
the matrix v{k) on the real axis according to 



(4.13) 
Set 



V = br hiT = 



-A^{k)R{k)e~^''s{k) 
1 



1 

A-2(fc)i?(fc)e2"9(fe) 1 







m6y\ 


ken'^ 


(4.14) 


m(i) = < 




ken^ 






, m, 


else. 



Note, that this deformation respects our symmetry condition (j2.19p . Evidently, 
the matrices 6l and bu have jumps on Eq. The new jump matrices v^^\k), that 
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yU 


' IC 






-—"^rr 


1 ia 




-• — ir 




► — nL 






yL 


(—ia 
i — ic 





Figure 4. The first deformation step 



correspond to m^^\k) on tfiis contour, are 

git(ff+-9-) 



(4.15) w(i)(/c)=<( 



{R- -R+ + x)A^'^e*^9++g-) Q-itig+-g-) J ' 



c'*(9+ -f- ) -R+ + x) A2e*(9++9- ) 

e-it(9+-ff-) 



Again Lemma 13.21 sliows that the off-diagonal entries vanish. 
Now set = Sf, \ Sq, that is 

= E;^, U S^, = [ic, ia] U [-ia, -ic]. 



fee E^, 



After the deformation (|4.14p the jump along the real axis disappears. Taking into 
account property (c) of Lcmma l4.1l we obtain a new Riemann-Hilbert problem 

(4.16) 
where 



TO^^^(A:) = mL^'(fc)w(i'(A:), m^^' (fc) ^ (1, 1), k ->■ oo, 



(4.17) t;(i)(fc) = <' 



' v{k), keTfuTf, 



A-^{k)x{k) e-"(9+-9-) 
git(ff+-g-) A2(fc)x(fc)\ 

g-it(ff+-S-) ) ' 

e'*^ j ' 



j = l,...,iV, 



1 

A-2(fc)i?(fc)e2«9(fe) 1 



1 -A2(fc)i?(fc)e-2its(fc) 
1 



Note that Img(fc) > (resp., lmg{k) < 0) on the contour \ {ia} (resp., 
\ {— ia}), and the corresponding matrices are exponentially close to the identity 
matrix except for small vicinities of the points ±ia. 
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Our next step of conjugation deals with a factorization of the jump matrices on 
the set Sac- To this end consider an auxihary scalar Riemann-Hilbert problem (cf. 
[53]): Find a function F{k) = F{k, ^) analytic in the domain C \ Sc and a constant 
/i(^) such that the following properties hold 



. F+ik)F^ik) = |x(fc)| for fc e S^, = [ic,ia] 



. F+ik)F^ik) - \xik)\-^ for e S^, - [-ia,-ic], 

• F+{k) = F_{k)h for A: G Sa = [ia, -ia], 

• F{k) ^ 1 as /c ^ oo and F{-k) = F~^{k) for fc e C \ E^- 

Note, that the last property allows us to use the function F as an entry of the 
diagonal matrix for a conjugation step. 

We construct the function F using the Plemelj formulas. In the domain C\ Eac, 
introduce the function 



(4.18) wik) = v/(fc2 + c2)(fc2 + a2)^ y^(o) > 0, 
and for fc e Sc set p{k) :— w(fc)-|_ — w{k)r. Then 

(4.19) p(-fc) = -p{k) for k G Eqc, p{k) = p(-fc), as A; e Sa. 
Set also 

log|x(fc)| 



(4.20) /(fc) 



p{k) 



Taking logarithms of the jump conditions and dividing them by p{k) we get 



^ 27ri yji^ s-k s-fc w(s)(s - A;)y J ' 

Properties (j2.14p and (|4.19p imply F{—k) = F^^{k). From this property, decom- 
posing the function in exponent with respect to k at infinity we conclude that 

(4^22) F{k) = 1 + + O (1 

,«3) „,o.±l-2r;!M..+iAr"'.'* 



27r w+{s) Ji^ w{s) 



where w{k) — ^ ik"^ + c^){k'^ + a?), 



fie iog|x(s)l ^„ 



(4.24) A(0 = 2i e M, and MO - e'^«) 

J— ia u;(s) 

Now we are ready to perform the next deformation-conjugation step. Introducing 
F{k) — F{k)A~^{k) we get again a function satisfying F{—k) — F^^{k) and hence 
the symmetry conditions of Lemma |3. II Moreover, observe x/i = \x\ on x/i = 
— Ixl on E^j, by (|2.14p . Using also condition (b) of Lemma 14 . 1 1 one can check, that 
on the contour E^^ the jump matrix w^^^ {k) can be factorized as 
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/ 1 



-i 
1 / V-i 



0^ 



1 ""2,+ 



where 
(4.25) 



D2{k) = 



F-^{k) 



F{k) 
From (|2.13p we conclude, that 





(7 



D2,-Gtik) 



-i 



L 



where 
(4.26) 



G^(yfc) 



TTi 
1 



G^(fc) 



1 



TTiF2A2 



Note, that the function TTi = ^j/^^'-.^ ^ has no jump along the real axis and is an 
even function in the domain C \ Ec- Introduce the symmetric domains 51^ and fl^ 
as depicted in Figure [5j Their boundary contours are oriented top-down. 





• — "V^rf 




-'I 


-~- — - — ir 


Ea 


CL 




E^ 




Q—ia 

' 5^ 

'—ic A 


" E^ 



Figure 5. The second deformation step 
Introduce the new function 
(4.27) TO(2)(fc) = m(i)(fc)G(fc), G(fc) = 



D2{k)G'^{k), kenY, 
D2(k)G'^{k), k € n^, 
D2{k), else. 



Applying again Lemma [3. II we arrive at a new RH problem: 



(4.28) 



^_^\k) ^m^^\k)v^^\k), m^^^(fc) -> (1,1), k^ 



(2)^ 
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where m^^\—k) = m'-^^k) 



and 



(4.29) 



/^Hk) = { 




u 



fc e E 



k e 


Ea , 


fee 




fee 




k e 





Here B = B + f (recall that F+Fr^ = e''^ on Ea, where A is defined by gS!])). 



Note that due to 



, (|2.13l) and the definition of F we get 
lim ^ and lim G^^ ^ 0, 

fc— >-ic fc— 7- — ic 



and, therefore 

(4.30) m(^)(fc) = (ci,C2)(A:Tic)"^/'*(l + 0(l)), as fc -> ±ic, where ci 7^ 0, C2 ^ 0. 

Now observe that on the contour 

E Uj^i(Tj^ U Tj^) U Ef U Ef U E^ U E^, 

all jumps are exponentially close to the identity except for small vicinities of the 
points ±ia as i — )■ cxd. To remove those parts one needs to solve the corresponding 
RH problem corresponding to the jumps on E^ U E*^ restricted to a small neigh- 
borhood of +ia (as well Ef U E^ restricted to a small neighborhood of — ia which 
however follows from the first by symmetry). Following the arguments from 
one can show that the contributions of these neighborhoods are negligible, that is, 

q{x,t) = qo{x,t) +0(1), 

where qo{x,t) is obtained from a "model" RH problem, where all jumps on E are 
discarded. We will solve this model RH problem in the next section. 



5. Solution of the model RH problem 

Consider the two-sheeted Riemann surface X associated with the function w{k), 
defined by (|4.18p . where we choose the standard branch of y/^ with the cut along 
the negative axis. The sheets of X are glued along the cuts [ic, ia] and [—ia, — ic]. 
Points on this surface are denoted hy p — {k, ±). The canonical homology basis of 
cycles {a, b} is chosen as follows: The a-cycle surrounds the points — ia, ia starting 
on the upper sheet from the left side of the cut [ic, ia] and continues on the upper 
sheet to the left part of [— ia, — ic] and returns after changing sheets. The cycle b 
surrounds the points ia, ic counterclockwise on the upper sheet. Moreover, consider 
the normalized holomorphic differential 



(5.1) 



duj = 27ri 



dk 



dk 



w{k) \7a w{k) 
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then duj = 2m, r ~ r(^) — j^dut < 0. Let 



J^expji 



be the theta function and recall that 9 is an even function, 0{—z) = 9{z), satisfying 
e{z + 27rin + T{i)l) = 6'(z) exp |-1t(0£2 „ 

Furthermore, let A{p) = J^^du) be the Abel map on X. Note that on the upper 
sheet, where p = (A:, +), it has the following properties: 

• A+{p) = -A^{p){ mod 27ri) for p e Ysac, 

• A+{p) - A_(p) = -r as p e 

• A{-p) = -A{p) + 7ri( mod 27ri) as fc G C \ Sc, p = {k, +); 

• A(ia) = — 1-( mod r), ia) = — ^ — 7ri( mod r, mod 27ri); 
. A((oo,+)) = f . 

Finally, denote by X = ^ + 7ri the Riemann constant associated with X . 

Identifying the upper sheet of X with the complex plane we introduce two func- 
tions 

(5.2) a\k) = 6 (A{k) + ^--K-'^y (^A(fc) + Z-+m-K-'-^ 

(5.3) f3\k) = e (^-A{k) + l-K-'-^y (^-Aik) + I+,ri-if-^ 

where 6 e K. will be determined later and A{k) — A{{k,+)) for k E C. 

Evidently, both functions and have zeros of order one (on X) at the points 
±ia. Moreover, 



/ 7ri 




/ 7ri 


if6\ 











(5.4) 

Due to the first three properties of the Abel map we get 

(5.5) a^ik) = ptik) and /3^(fc) = a^_(fc) for ke^ac^ (S^, U S^J. 

(5.7) a^-k) = /3^(fc) for fc e C \ S^. 

Now introduce the function 



(5.8) 7(fc) = 

defined uniquely on the set C \ Sc by the condition arg7(+0) = 0. This function 
satisfy the jump conditions 

..q^ 7+(fc)=i7-(fc), fcGS^c 

^ ' 7+(fc) = -i7-(fc), fcGS^e- 

Combining (|5.4p - (|5.9p we conclude that the vector 

(5.10) ^(3)(fc) . (,ik)^^^m^, ,ik)^^m^] 

V aO(fc)aS(c3o) /?0(fc)/3S(oo)y 
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solves our model problem 



(5.11) 
where 



(5.12) 



^^^\k) ^m^^\k)v^^'>{k), m^^)(fc) -> (1,1), k 



(3)1 




The symmetry condition 
(5.13) m(3)(-fc) 



1 

1 



k G Saj 



m^^Hk), keC\^c 



is also fulfilled due to (|5.7p . 

Moreover, both components of the vector- valued function m'^'^\k) are bounded 
everywhere except for small vicinities of the points ±ia, ±ic, where they have 
singularities of the type (fc — C)~^^^; C € {ic, ia, — ic, — ia}. 

In summary we get 



(5.14) mf\k) 





e 


{m - i 






e{A{k)-m)9{A{k))e 


[2 2 ) 


^( 





(5.15) (fc)= ^ 



-A(fc) - ) 02 (-i) 



(-A(fc) - iTT) (-A(fc)) - if ) (ii + if) 



We are interested in the terms of order as fc — > +ioo. To this end let 



(5.16) 



ds 



V/(C2-S2)(a2_s2) 



< 



be the normalizing constant from the Abel integral. Since w{k) — fc2(l + o(l)) as 
k — > +ioo, we infer 



A{k) - A{+oo) = A{k) 



k 



O 



and 



Proceeding in the same way with the other theta functions and taking into account 
that 7(fc) = 1 + 0(fc^2) for large k, we get 
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where 

E{u) = E{u, = log {e{u)e{u - m)) . 

To check that these asymptotics are real- valued on the imaginary axis, recall that 
^(") = ^3 (2^), where (cf. [2]) 

= I ^1) = 1] exp{(mVi + 2mv)Tii), n = ri(0 = ^ G iM+. 

ZTTl 

Then i9(«) = ^,^e,{v) and 

where 



(5.17) Eiv) = E{v,0 = ^|;log (^03(«)e3(^' - I) 

is a real-valued function for w G M. 

Since outside small vicinities of the points ±ic, ±ia the solution of the model 
problem tvP'^k) approximates the solution ra'^'^^k) of the problem (I4.28l) - (|4.29p it 
remains to trace back our deformation and conjugation steps: 

"^y^iy ^{k))\ d{k,t))\ A-i(A:)i^(fc)y 

where the asymptotic behavior of d(k,t) and F{k) are given by ()4.8p and (|4.22p - 

(|4.24l) . Therefore, with i|j = ^, we have 

(5.18) 

.(., . (1, 1) + ("(^)-^||-^) + ^lil±i(i) ) (1, -1) + o (i. 

This formula together with (j2.16p gives us the asymptotic behavior of q{s, t)ds. 
To get the asymptotic behavior of q{x, t) we differentiate (|5.18p with respect to x, 
taking into account that for any smooth function p{£_) one has j^p{£,) = 0(1). 
Thus we obtain 

where 

(5.20) E'{v) = -f (0^ log (^3 {v I ri(C)) 9, (^{v - 1) | n(0 
and 

(5.21) f(e) = i(^|"^''\(c^-.s^)(a^^ . 

Formula (|5.19p can be simplified using the following formula for summing theta- 
functions (14 formula (1.4.3)) 



73 



{z + w\ -)93iz ~w\-)^ 03(2z I t)03(2u; | t) + ^?2(2z | r)92i2w \ r), 
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where 

62{z I t) = ^ exp{7ri(TO + -)^r + 27ri(m + -)z}. 

Since ^2(5 I ''") = 0, we see 

0,{u\^^)9^{u-\C-) = e,{2u-\\r)e,{\\T) 
and the last formula implies that 

log f^a {v I ri(0) 03 f - ]:) \ ti(o')') = \ogd:,{2v - ]- \ 2ti(0) + /(O- 



Substituting w = i — and taking into account (|5.19p as well as the 

estimate = 0{\) we get 



Theorem 5.1. Assume (|1.3l) - (|1.4p . T/ien m i/ie domain ~Qcr + e < f < 4c^ — £1 
t/ie following asymptotical formula is valid 

(5.22) + 

-ffere 03(w) = Qj,{v \ t(^)) wif/i r(^) = j\^duj, where duo is the normalized holo- 
morphic differential (|5.ip . T/ie function T{^) > is defined by ()5.2ip bkc? 



^(0 



2 

iog|(T(is)ri(is)| ^ / r(«) 



Ja(n J(c^-s^)(s^-a(m \J- 



are real-valued functions. In all formulas for T, B, z and A the positive values of 
the square roots are taken. 

6. ASYMPTOTICS IN THE DOMAIN X < —Qc'^t. 

To study the asymptotical behavior of q{x, t) in the domain we use the RH 
problem, associated with the left half axis. Namely, we consider the spectral data 
and the Jost solutions as the functions of the parameter ki = \jW^-V<? . Then 
the continuous spectrum of the operator i?(0) coincides with the set Im(A;i) ~ 



and the discrete spectrum is located at the points i^ij — ^ (recall that 

k| > (?\ Introduce the vector- valued function 
(6.1) 

r (Ti(fci,i)0(fci,a;,i)e-"=^^ 0i(/ci, x, i)e"^i-) , fci e C^, 
^ ' ((/.i(-fci,x,i)e-"=i- ri(-fci,t)0(-fci,a;,t)e"=^-), fci e C^, 

where C'^ := {fci : Im(fci) > 0}, {kx ■ Im(/ci) < 0}. 
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This function has the fohowing asymptotical behavior 

1 



(6.2) m{ki,x,t) ^ + 



2ifci 



{qiy,t)+c')dy (1,-1) + 



Theorem 6.1. Let {Ri{ki), ki e M; (ki.j, 7i.j), 1 < J < N} be the left scattering 
data of the operator H{0). Let (resp., T^') be circles with centers in iKi.j 
(resp., —iKij) and radiuses < e < imin^j^jKij — ki^ = 0. Then 

m(fci) = m{ki,x,t) defined in ()6.ip is a solution of the following vector Riemann- 
Hilbert problem. 

Find a function m(fci) which is holomorphic away from the contour U^]^(Ty U 
T^) U R and satisfies: 

(i) The jump condition m+(fci) = m_(fci)w(fci) 



l-|i?i(fci)|2 
i?i(fci)e**i('=i) 



-i?i(fci)e 
1 





i?i(fci)e 



t*i(fei) 



-i?i(fci)e-**i('=i: 

1 



(6.3) w(fci) < 



1 



-1,3)- 



1 



fci e R \ [-C, c] 
fci e [-C, c] 



1 

1 



(ii) the symmetry condition m{~ki) — m{ki 

(Hi) the normalization condition ImiK^oo m{iK) — (1 1) 
Here the phase $i(fc) — $i(fci,x,t) is given by 

-8ikl + 12ic^ki~24i^ki, ^ 



(6.4) 



$i(fci) 



X 



Proof. The proof of this theorem is similar to the proofs of Theorem 
131 It is based on formulas (EH]), (EJT]), (I^TU)) . the formula 
4>(ki, x,t) eM. for ki E [— c, c] and the relations (cf. [15] ) 



and Lemma 
-fci,x,i) = 



i?i (fci , i) = i?i (fci , o)e-4it(-c^+2fe^)/ci ^ ^2 . ^ ^2 . 



3k - Ki jt — Ac Ki 



□ 



Denote by i/ci^o = ± 



^ the stationary phase points of $1, that is, the 



zeros of the equation <^'i{ki) = 0. In the present domain ^ < — ^ we have ki Q > c 
and the signature table for Re$i is shown in Figure [5] 

First of all we observe, that the jump matrices corresponding to the discrete 
spectrum are exponentially close to the identity matrices as t — >■ 00. Therefore, 
unlike in the previous cases we do not need a conjugation step for them. Moreover, 
since the parameter k will not appear in the remainder of this section, we will write 
k in place of ki to simplify notations. 
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Figure 6. Sign of Re($i(fci)) 
From (Pl^ it follows that 1 - |-Ri(fc)p = for fc e [-c, c] and hence 



Now following the usual procedure [12], [TS] we let d{k) be an analytic function in 
the domain C \ (M \ [— fci.o, fci.o]) satisfying 

= d_(fc)(l - for fc e M\ [-fci,o,fci,o] andd(fc)^l, k^oo. 

Then by the Plemelj formulas 

For the smooth steplike initial data q{x, 0) G C"(R) the reflection coefficient satisfies 
Ri{k) — 0{k"'^^) (cf. 15 ). Moreover, in the domain M\ [— fci.o, fci.o] we have 
|i?i(fc)| < 1. Therefore, the integral under the exponent is well defined. Since the 
domain of integration here is even and the function log(l — |i?ip) is also even, we 
obtain fc) = d^^(fc) and the matrix 

^d-i(fc) 
rf(fc)^ 

satisfies the symmetry conditions of Lcmma l3.1l Now set rh{k) = m{k)D{k) and the 
new RH problem will read rhj^(k) = rri^{k)v{k), where m(fc) — )• (1, 1) as fc — > oo, 
'm{~k) = rh{k) ( ? q ) and 

' A^_{k)Al{k), fceR\[-fci,o,fci.o] 

(6.7) i{k)=\ B^{k)B^{k), fee [-fci,o,fci,o] 

_ D-^{k)v{k)D{k), fc e Uj(T^^UTj^), 

where 

(6.8) A^{k) = { J, kenfyjn^, 

V(i-ifli(fc)P)d2(fc) V 



(6.6) D{k) 



/-, d"(fc)fli(fc)c-"^i"') \ 

(6.9) A^{k)=\^ (i-l-Ri(fe)l^) j, fcef^Puf}^, 



(6.10) ; -rf^(W(^)e-^^'=^V fc,^. 
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(6.11) 



ken'-. 



Here the domains fif, nf , f2^, f2^, fi^, fi^, are bounded by the contours Sf, 
EY, S^, E^, S^, E^, as shown in Fi gure [71 All contours are oriented from left to 



Ro($i)>0' 



Rc($i)<a- 



Ro($i)<0 




Rc(*i)>0 



Figure 7. Contour deformation in the dispersive region 



right. They are chosen to respect the symmetry fc i— > — fc and are inside the strip 

L inside the domain, where 
Ri{—k) in these domains. 



I Imfc| < min{^, ^^j^} below the discrete spectrum and inside the domain, where 



3.12) 



i{k) = < 



kenyu n^, 
k e nj^, 
k e n^, 

else. 



Ri{k) has an analytic continuation. We also set Ri{k) 
Now redefine rh{k) according to 

rh{k)A^{k), 
rh{k)A^{k)-^ 
m{k)B^lk), 
m{k)B^{k)-^ 

77l(fc), 

Now the function 77i(fc) has no jump asfc g K and all evidently defined jumps on 
contours Ef, E;^, E^, E^, E^, E^, u|Li(Ty UTj^) are exponentially small with 
respect to t outside of small vicinities of the stationary phase points fci^o and — fci^o- 
Thus, the model problem here has the trivial solution rh{k) = (1,1). For large 
imaginary k with |fc| > ki^at + 1 we have fh{k) — m(k) and consequently 

m(fc) = m{k)D-^{k) = {d{k),d-\k)) 

for sufficiently large k. By (|6.5I) 



rf(fci) = l + -^(-- 

2lfcl \ TT 



iog(i-ii?i(s)nds 



Comparing this formula with formula (|6.2p we now can derive the asymptotics using 
Theorem A.l from [2^ following literally the argument in Section 5 of [TS]: 

Theorem 6.2. Assume p.3p and (jl.4p . Then the asymptotics in the similarity 

region, j + 6c^ < — C for some C > 0, are given by 

(6.13) 



q{x,t) 



4i^(fci^o)fci,o 
it 



sin(16i/c?_o " '^(^i.o) log(192ifc3 o) + (5(fci,o)) + 0(i"") 
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2t '^"'^ 



for any 1/2 < a < 1. Here ki^ — \J ^ ~ j 
(6.14) i.(fci,o) = - log(l - |i?i(fci.o)n, 

ZTT 

(^Ifci.o) = J - arg(i?i(fci,o)) + arg(r(iz/(fci^o))) 

1 /■ i-|i?i(C)P A 1 



TT 



(6-15) - - / log , :^^rfc 



fei,o,fei,o] Vl ~ l^l(fcl,o)IV 



Appendix A. Inverse scattering transform on steplike backgrounds 

The purpose of this Appendix is to prove some facts from scattering theory used 
in this paper. Most of the properties hsted here are vahd for a much wider class 
of potentials then those satisfying (|1.3p . namely, for continuous potentials with a 
finite second moment: 

(A.l) / {I + x^){\q{x,t)\ + \q{-x,t) + c^\)dx <oo. 

Jo 

We start with 

Proof of Lemma \2.S[ We will omit the dependence on t for notational simplicity. 
Let (j>{k,x) and (f>i{k,x) be the Jost solutions (|2.4I) of equation (|2.2[) . normalized 
by (|2.3p . According to (|2.6p the right transmission coefficient T(fc) is defined by 
formula (|2.10p . Our first step is to compute its asymptotics as /c cx) up to a term 
o (-j:). Since the Wronskian (|2.7p does not depend on x, we evaluate it at x = 0. 
Under condition (lA.ll) the integrals in (|2.4p can be integrated by parts one time 
and then differentiated with respect to x. We get 

0l(fc, 0) - 1 - ^^^^T^ + ° ' "^'l^^' = + + 

(^(fc, 0) - 1 - + o ; 0'(fc, 0, t) = ik- K{0, 0) + o(l). 

Since ^ = l + C)(p-)we further infer 

k fci 

iy(fc, i) = ifc - — Ji:i(o, 0, t) - ii:(o, 0) + ifci - -t-k{o, o) 

fci fc 
- Xi(0, 0) + o(l) = 2ifc - 2{K + A'i)(0, 0) + o(l). 

Thus, 

(A.2) T(fc) = l + ^iM)±^im+,fl 

ifc' \k 

Next, since fci = ^ + |^ + O ( Jj-) , 
(A.3) 

\ ifci \fc 

=e-%^(^+°(^))^ f 1 - ^^^) + o f-) = fl + — 
\ ifc / \k J \ 2ifc 

xU-^lii^UoaUl + fe^^i^ + o 



X 



\k ) \k J ifc- V fc 
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From formulas (12.51) it follows, that 



d 

— {Ki{x,x)+K{x,x)) = — 

and therefore 

2 

(A.4) Ki{x,x)+Kix,x) = ^+Ki{0,0)+K{0,0). 

Combining (|A.2|) - (|A.4p we get 

kx , K{x,x) ( 1 



T(A;)0i(fc, x)e"=^ - 1 + + o[- 



-ikx _ 1 K{x,x) ^ f 1 



On the other side, 

0(fc,a;)e . ' " \ k 

which finishes the proof. □ 

Proof of Theorem \2.4\ We begin by checking that the jump condition for the vector 
m(A;, x, t) defined in (j2.15p has the form (j2.17p . Since our further considerations are 
mostly algebraically, we omit the variables x, t and sometimes also k in notations 
whenever possible. 

Consider fc g S = R. Let a ) be the unknown jump matrix. Since T{~k) = 

T{k), 01 (—fc) = <pi{k) for fc e S, the entries of m satisfy 

1 (pie = (pe a + 1 (pie 7, (pe = (pe p + 1 (pie d. 

Multiply the first equality by e~''^^, the second one by e''^^, and then conjugate 
both of them. Abbreviating 

(A.5) 7' = 7e-2"=^ and^' = ^e2"=^, 

we finally get 

a(l> — T(/)i — T^'(j)i, T5(j)i = — 
Now divide the first by T and compare both with (|2.6p . This shows 5 = 1, — /?' = 
R{k,t), a = TiT, and -Y= = Ri. Applying ((X5|) . (PTT|) and the evolution 
formula for R{k,t) from Lemma [2.11 5., we get v{k) for fc g E. 

Now let k G C C"*" implying that we have to work with the upper case in 
(|2.15p . Since the function (j){k) is real- valued and has no jump on this set, the 
equations for the entries of the jump matrix read as follows (assuming k G [0, ic]-^) 

T^j, ^ikx rri 1 Akx , i^ — ikx A^—ikx rji I Akx n I i^ — '\kx 

1 (pie = 1 (pie a + (pe 7, (pe 7 = J (pie p + (pe 

From the last equality f3 = 0, 5 = 1. Abbreviating 7e~^''^^ =: 7' and divide this 
equality by —T. By virtue of (|2.12p we get 

— -7' 
-I- (/)ia = -^0- 

Comparing this equality with the first of the scattering relations (12. 6p we obtain 
a = 1 and 7 = ~Ti{k,t)T{k,t)e^'^'' . Since -TiT = -^\Ti\^ for k G S^, the 
corresponding formula from item 5. of Lemma l2 . II establishes the formula for w(fc), 
k&YH. □ 
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